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SUMMARY 


The  following  analytical  study  examines  the  application  of  digital  technology 
to  seismometer  measurements  in  borehole  installations.  Current  practice  is  to 
generate  within  the  seismometer  system  an  analog  output  projtortional  to  seismic 
input  and  transmit  this  signal  in  analog  form  to  the  top  of  the  borehole.  1 he  inherent 
disadvantages  of  transmitting  low  level  analog  signals  over  long  transmission  lines 
underline  the  importance  of  converting  the  signal  output  to  digital  format  the  bottom 
of  the  borehole. 

Several  methods  of  digitizing  the  seismometer  output  were  examined.  Two 
methods  involved  closing  a digital  torque-to-balance  loop  around  the  seismometer. 
Both  ternary  (three  state)  and  binary  (two  state)  loop  implementations  were  analyzed 
and  simulated  using  an  IBM- 360  computer.  The  ternary  loop  was  found  unacceptable 
since  it  resulted  in  the  generation  of  an  input  dependent  noise  signal  riding  on  the 
desired  output.  This  noise  could  not  be  filtered  since  it  could  vary  over  a wide 
frequency  band. 


The  binary  pulse  width  modulated  loop  performed  well  in  the  computer 
simulation.  However,  use  of  the  binary  loop  resulted  in  the  generation  of  a high 
ripple  signal  summed  with  the  desired  error  signal.  The  presence  of  this  signal 
put  impractical  requirements  on  the  electronics  within  the  control  loop,  thus 
eliminating  the  binary  loop  from  further  consideration. 

The  two  remaining  implementations  studied  involved  the  use  of  an  analog 
restoring  loop  in  conjunction  with  an  analog-to-digital  (A/D)  converter.  The  A/l) 
converter  was  assumed  to  be  either  inside  the  control  loop  or  external  to  the  'oop. 
The  v anfiguration  with  the  A/D  conversion  inside  of  the  loop  proved  far  superior 
and  resulted  in  an  overall  system  performance  which  was  quite  good. 

The  worst  case  noise  caused  by  the  A / D converter  was  only  slightly  higher  than 
the  noise  in  the  current  analog  loop  configurations  (llpV  rms  (0-10  Hz)  vs  4 fN  rms 

(.02-1  Hz)). 


The  A/D  converter  provides  an  lb  bit  dynamic  range  and  retains  the  simplicity 
of  the  add-on  external  approach.  The  electronics  design,  discussed  in  detail  in 
Chapter  5,  can  be  built  at  reasonable  cost  and  requires  no  overly  complex  circuitry. 
Errors  due  to  nonlinearities,  temperature  variations  and  other  instabilities  would 
be  reduced  by  the  loop  gain  and  could  be  made  negligible.  In  the  case  of  the  outboard 
A/D  converter  there  would  be  no  reduction  of  these  errors  and  additional  effort 


vii 


would  be  required  to  adequately  compensate  for  them. 


The  study  concludes  that  the  inboard  A/ 1)  converter  provides  an  attractive 
alternat've  and  recommends  that  an  18  bit  A/I)  converter  be  interfaced  vutha  Geotech 
36000  instrument  and  the  resultant  performance  be  evaluated  ata  "seismicallyquiet" 
site. 
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ClIAt'di  It  1 
I \ Tift  1 1 )l  C 1'ION 

hi  recent  years  the*  developn  ent  of  sophisticate  i digital  control  and  encoding 
techniques  lias  rr  ade  availabh  to  relate  1 iiscipl  ines  the  numerous  advantages  inherent 
in  such  techniques.  At  Draper  Laboratory  digital  control  loo|  s have  heen  put  to 
large  scale  use  with  a variety  of  inertial  grade  ln^trun  ents.  1 or  oxan  pie,  a ternary 
force-rcdialance  loop  is  used  m the  1‘oseidon  guidance  svsteit  with  the  Draper 
Laboratory  size  1 (i  I'M  I’ll.  I'hc  loop  provides  quantized  increments  of  velocity 
data,  ‘.linnlar  work  using  binary  and  pulse  width  modulated  control  loops  with  gyros 
and  accelerometers  has  been  done  m hardware  lesign  optimization  studies,  i he 
effort  in  these  areas  has  been  to  achieve  a high  level  ol  resol  it  ion  witii  extren  oly 
good  long-tern  pulse  weight  stability. 

1 he  wo  rk  described  in  the  present  r«qort  is  an  alien  | t to  apply  ligital  technology 
to  seismometer  svsterns.  The  candidate  ins’run  ent  used  m the  analysis  was  the 
(ieo'ech  afiOOU  horizontal  sensing  seismon  eter.  I'his  instrun  ent  w is  n.odelled  with 
both  ternary  and  binary  pulse  width  modulated  control  loops,  in  i idition,  hi  analog 
configuration  was  evaluated  employing  either  mlioard  (within  the  feedback  loop)  or 
outboard  A/D  conversion  of  the  analog  out)  it.  Compiler  simulations  were  used  to 
verify  linearity  and  loop  stability  md  to  exan  ine  tne  loo|  response  to  sin  ulated 
inputs. 

Chapter  2 contains  a brief  description  of  the  (ieotech  MtiOOfl  instrument.  A 
basic  analytical  model  is  derived  and  a brief  sun  nun  of  instrument  operating 
specifications  is  presented. 

Chapter  3 presents  a theoretical  analysis  of  the  various  digital  implementations 
used  with  the  basic  seismometer  model.  Most  of  the  chapter  is  devoted  to  a detailed 
design  of  the  pulse  width  modulated  binary  configuration,  since  this  configuration 
appeared  to  lie  the  onl\  viable  fullv  digital  in  plen  entation,  A prime  disadvantage 
of  this  configuration  is  the  large  ripple  component  riuing  on  the  desired  dc  error 
signal.  This  ripple  was  not  at  first  considered  to  be  a serious  problem;  however, 
when  the  electronics  design  was  done  in  detail  (as  shown  in  Chapter  ;>)  it  was 
discovered  that  the  presence  of  the  large  ripple  component  necessitates  use  of  a 
sample  and  hold  circuit  with  an  accompanying  large  phase  lag.  I lie  total  loop 
configuration  becomes  unstable  and  extremely  difficult  to  implement,  if  at  all  |>ossible. 

The  latter  iiortions  of  Chapter  3 examine  an  analog  configuration  with  A D 
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conversion  within  the  control  loop.  Though  this  approach  was  the  last  chosen  for 
analysis,  it  exhibits  definite  advantages  over  the  other  approaches  and  an  attractive 
option  for  further  development. 


Chapter  4 details  a basic  approach  to  the  electronics  design  taking  into  account 
the  seismometer  design  goals,  while  attempting  to  maintain  a final  design  concept 
which  remains  both  feasible  and  economical.  Hie  detailed  electronics  design  is 
presented  in  Chapter  n.  The  approach  taken  is  general  in  form,  applying  in  varying 
degrees  to  all  of  the  digital  implementations  described  in  Chapter  3.  The  result  of 
the  studx  is  the  Highlighting  cf  some  of  the  problems  associated  with  the  various 
implementations  and  a selection  of  the  inboard  \/l)  conversion  as  that  approach 
appearing  to  offer  the  most  promise.  \n  analysis  of  the  noise  resulting  from  the 
A/P  converter  is  presented  in  Chapter  fi. 


Chapter  7 presents  conclusions  resulting  from  the  study  and  recommendations 
for  future  effort. 
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CHAPTER  2 

DESCRIPTION  OF  THE  36000  SEISMOMETER  OPERATION 
2.1  Physical  Description  and  Operation 

The  36000  seismometer  is  a horizontal  sensing  pendulous  type  device  capable 
of  measuring  long  period  seismic  disturbances  (periods  of  10  to  100  seconds).  The 
seismometer  assembly  consists  basically  of  a mechanical  spring  restrained  pendulous 
mass  attached  to  the  case  or  frame  of  the  seismometer  and  free  to  pivot  in  a single 
plane.  Provisions  are  available  for  both  dissipative  damping  (internal  friction  and 
external  electronics)  and  nondissipative  damping.  In  addition,  an  electrodynamic 
transducer  (moving  coil)  is  used  for  electronic  feedback  spring  control. 

Motion  of  the  pendulum  or  mass  with  respect  to  the  case  is  measured  by  a 
displacement  transducer  consisting  of  a capacitance  bridge  arranged  such  that  the 
capacitance  varies  linearly  with  motion  of  the  mass.  The  bridge  is  excited  by  a 
carrier  signal  of  constant  magnitude  and  frequency.  Movement  of  the  pendulum  is, 
therefore,  translated  into  an  amplitude  modulated  signal.  A photograph  of  the  36000 
assembly  is  shown  in  Figure  2.1-1. 
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Fig.  2.1-1  Geotech  36000  Seismometer 


A block  diagram  of  the  system  is  depicted  in  Fig.  2.1-2.  The  carrier  signal 
is  amplified  by  an  IF  amplifier  and  then  demodulated  to  yield  an  analog  signal. 
This  signal  is  measured  and  fed  back  to  the  electrodynamic  force  generator  to  provide 
restraint  on  the  motion  of  the  pendulous  mass.  Additional  electronics  divide  the 
signal  into  low  and  high  frequency  intervals  (less  than  and  greater  than  0.02  Hz) 
and  provide  necessary  loop  compensation. 


HIGH 

SEISMIC'  FREQUENCY 

DISTURBANCE  SIGNAL 


LOW 

FREQl KNCY 
SIGNAL 


Fig.  2. 1-2  Basic  Block  Diagram  of  Seismometer 


2.2  Modeling  of  the  Seismometer 

The  suspension  block  depicted  in  Fig.  2.1-2  was  modelled  as  a garden  gate 
type  of  device  sensing  accelerations  in  the  horizontal  plane  with  the  axis  of  the 
pivot  vertically  oriented  (see  Figure  2.2-1). 

Let: 


M = mass  of  the  boom  (kg) 
k = radius  of  gyration  (m) 

0 = angular  displacement  of  the  mass  (rad) 

d = damping  constant  (N-m-s/rad) 

K = spring  constant  (N-m/rad) 
r = pivot-center  of  gravity  separation  (m) 
x = horizontal  acceleration  of  the  case  (m/s  ) 
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Summing  of  the  torques  about  the  pivot  yields. 


■)..  „ r>-i) 

Mk  e + de  + KG  - Mrx  =0  v" 

Taking  LaPlace  transforms  of  both  sides,  the  transfer  function  relating  angular 
displacement  to  input  acceleration  is  obtained. 


0(s) 

X (s) 


Mr/K 

Mk2s  /K  + ds/K  + 1 


(2-2) 


The  capacitance  bridge,  IF  amplifier,  demodulator  and  other  amplification 
were  modelled  as  CqKq. 

Dynamics  of  the  associated  electronics  which  divide  the  analog  signal  into 
low  and  high  frequency  bands  yield  the  following  transfer  function: 

(TjS/U  + Kq)  + U/Ciy  + 1)  (2-3) 

where  T^/1+Kq  corresponds  to  50  seconds  and  to  10,000  seconds. 
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The  compensation  block  provides  lead  phase  to  compensate  for  the  lag  phase 
introduced  by  the  suspension  and  the  signal  divider.  A resistor  and  capacitor  in 
parallel  give  the  requisite  lead  for  the  loop  and  current  levels  for  the  torquer. 
The  transfer  function  relating  current  to  voltage  is  given  by: 

i (s)/e (s)  = <T2s  + 1)/R  (2‘4) 


where: 

R = resistance 

and  T2  corresponds  to  about  1.8  seconds  and  compensates  for  the  180°  phase  shift 
introduced  by  the  suspension  at  its  natural  frequt  ncy  of  5 seconds. 

9 

Finally,  the  torquer  is  represented  as  a simple  gain  with  the  units  of  m/s  /A. 

This  linear  model  of  a seismometer  loop  provides  the  point  of  departure  for 
the  analysis  of  other  configurations  discussed  in  Chapter  3. 

2.3  Noise  Specification  and  Instrument  Parameters 


The  design  of  a good  seismometer  requires  that  the  earth  motion  equivalent 
of  thermal  agitation  noise  introduced  by  the  suspension  and  the  associated  electronics 
should  be  20  db  below  seismic  background  noise. 

At  a quiet  deep  mine,  suchasthat  located  at  Queen  Creek,  Arizona,  the  minimum 
background  noise  level  (at  0.035  Hz)  is  3.4  x 10'10  m/s2/  Hz.  The  maximum  noise 
level  introduced  by  the  seismometer  should  be  20  dB  down  from  this  value  or 
3.4  x lO"11  m/s2/  Hz.  This  is  considered  to  be  a design  optimum. 


Table  2.3-1  presents  a tabulation  of  parameters  and  specifications^  for  the 
36000  instrument.  Tbe  noise  level  of  the  Geotech  instrument  equals  4 x 10  m/s 
(0.02  to  1 Hz). 


The  maximum  long  period  output  signal  level  stated  in  the  Geotech  specification 
(0.02  to  1 Hz)  corresponds  to  lO-4  m/s2  input  to  the  seismometer.  If  the  minimum 
signal  to  be  detected  is  at  the  noise  level  of  the  instrument,  then  a dynamic  range 
of  2 x 10^  is  required.  This  value  is  equivalent  to  a 21  bit  digitization.  The 
specifications  further  call  for  lO^o  linearity  over  the  dynamic  range. 


/ 
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Table  2.3-1 


Parame ters  and  Specifications 


Parameters 

Inertial  Mass 
Radius  to  Center  of  Mass 
Spring  Constant 
Moment  of  Inertia 
Dissipative  Damping 

Specifications 

Input 

Equivalent  Noise  (0.01  to 
1 Hz) 

Maximum  Input  Levels 
Output 

Noise  Level  (0.02  to  1 Hz) 

Sensitivity 

Offset 

Maximum  Level 
Response 


for  the  36000  Seismometer 


0.  364  kg 
0.0432  m 

_ O 

1. 0 x 10  ' N-m / rad 
6.  33  x 10’4  kg-m2 

_ C. 

1.17x10  ' N-m-s/rad 


21  0 2 

1.  56  x 10  (m/s")  /Hz,  max. 

-2  2 
±1  X 10  mis 

1 . 6 x 10"11  V2/ Hz,  max. 

1 x 105  V/m/s2  ‘5% 

+ 5V,  max. 
ilOV 

Flat  to  acceleration  input  from 
0 to  1 Hz,  down  12  dR/octave 
above  1 Hz. 


V-  - < * 


CHAPTER  3 


ANALYSIS  OE  FOUR  DIGITAL  SEISMOMETER  CONFIGURATIONS 
3.1  Introduction 

I our  digital  seismometer  configurations  were  analyzed  with  respect  to  their 
ability  to  produce  a linear  response  in  an  environment  of  accelerations  ranging 
from  10" 10  to  10‘2  m/s2  with  peri  ds  of  13  to  43  seconds.  It  should  be  noted  that 
the  total  dynamic  range  of  the  Geotech  feedback  loop  is  10  . Each  of  the  outputs  in 
the  current  Geotech  configuration  (0-1  Hz;  0.02-1  Hz)  has  a dynamic  range  of  10  . 
The  initial  study  presented  here  examines  a pulse  torque  feedback  implementation 
giving  a 106  dynamic  range  while  utilizing  the  output  voltage  range  employed  by 

Geotech. 

The  point  of  departure  for  the  overall  analysis  of  digitized  seismometer 
systems  was  the  analog  configuration  presently  used  by  Geotech.  Necessary 
modifications  were  made  to  this  configuration  to  achieve  the  required  digital 

performance. 

The  linear  configuration  for  an  analog  accelerometer  system  was  derived  in 
Chapter  2 and  is  shown  as  the  block  diagram  in  Tig.  3.1-1. 


Fig.  3.  1-1  Linear  Model  of  Analog  Seismometer  System 


where  the  variables  are  again  defined  as  in  Chapter  2,  with  the  addition  of; 

CQ  * AC  amplifier  gain  (V ! rad) 

K.  = DC  amplifier  gain  (V/Y) 

T = time  constant  associated  with  DC  amplifier  (seconds) 
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S LaPlace  transform  operator 

1 } time  constant  associated  with  stabilization  (seconds! 

C'j  torquer  constant  (N/A) 

R resistance  (\  A) 

I he  system  of  big.  H.l-1  generates  continuous  information  relating  to  ac- 
celeration. I our  alternatives  which  provided  a digitized  measure  of  , he  displacement 
x were  studied 


M 2 I erna  rv  Nystem 


■pile  ternarv  svstem  is  essentially  a ratt  measuring  device.  Signal  derivative 
information  is  converted  to  a pulse  rate.  The  pulses  activate  the  torquer  of  the 
seis.non  eter  and  a count  of  the  pulses  yields  velocity  information.  I igure  3.2-1 
shows  the  operation  ol  a typical  ternary  system. 


l ie-  M . 2-1  Steady  Stat<  T«  rnary  ( tperation 


I he  following  discussion  pertains  to  signals  of  positive  slope,  \ sirnplt. 
extension  may  be  made  to  signals  of  negative  slope. 

Let  e be  a signal  with  positive  slope  v.  Let  d be  a threshold  activating  voltage. 
When  the  signal  level  e reaches  d.  a pulse  of  amplitude  -1  is  subtracted  from  slope 
v for  a time  intervals.  Thu  system  then  resumes  "pulseless"  operation  uu  U the 
next  threshold  is  reached.  This  operation  is  shown  in  big.  ^2-2.  Dotted  lines 
imply  logical  connectors. 
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Fig.  3.  2-2  Block  Diagram  of  Ternary  Configuration 

To  determine  the  relationship  between  pulse  counts  and  v assume  steady  state 
(constant  v).  Let  T be  the  time  between  the  beginning  of  the  last  pulse  and  the  time 
that  threshold  d is  reached. 

For  steady  state,  an  examination  of  Fig.  3.2-1  shows  that: 

V(T  - l)  = "(v  - IK  or  v = tt/T  (3-1) 

i.e.,  the  time  between  pulses  is  inversely  proportional  to  v.  As  the  pulse  frequency 
f is  equal  to  1 /T, 


A counting  of  the  pulses,  therefore,  yields  velocity  information. 

A major  drawback  to  ternary  implementation  maybe  inferred  from  Lq.  3-2. 
The  frequency  of  pulses  exciting  the  torquer  is  a function  of  the  slope  of  the  signal. 
As  the  mass  responds  essentially  to  the  dc  value  of  the  pulse  train,  the  pulse  tram 
less  the  dc  value  is  noise  corrupting  the  performance  of  the  system.  This  noise 
can  not  be  simply  filtered  as  its  frequency  is  derivative  dependent. 

As  an  example,  assume  that  v is  required  to  a resolution  of  1 part  in  250,000. 
To  avoid  noise  frequencies  near  the  natural  frequency  of  the  mass  (approximately 
5 seconds).  T should  be  less  than  one  second.  This  in  turn  places  a constraint  on 
the  pulse  width.  J.  to  be  less  than  4 us.  The  noise  frequency  would  then  range  over 
frequencies  from  1 Hz  to  250,000  Hz  as  a function  of  the  magnitude  of  v. 
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V. 


3.3  lUn.tr>  System 


| iU!  binary  system  consists  of  a scheme  wherein  the  tnrquer  ol  tile  ac- 
celerometer is  excited  by  alternating  pulses  of  fixed  magnitude  with  the  constraint 
that  the  sum  of  tile  periods  of  the  negative  pulse  (T())  and  the  following  positive 
pulse  be  a constant,  1'  (see  I ig.  3.3-1). 


Fig.  3.  3-1  llinary  System  Fulse  Pattern 


T 

“Hie  dc  feedback  signal  to  the  tr.rquer  is  then  proportional  to  j-  T0.  Linearity 
is  maintained  by  making  I. TQ  propirtional  to  some  activating  voltage.  A system 
incorporating  binary  feedback  is  shown  in  I ig.  3.3-2. 


Fig.  3.  3-2  Block  Diagram  for  a Binary  Seismometer  System 


A compensation  network  considerably  more  complicated  than  that  of  the  linear 
system  is  required  because  of  the  delay  inherent  in  the  binary  implementation.  In 
order  to  effect  full  range  operation  (i.e.  a positive  pulse  lasting  for  a period  1 ),  it 
is  necessary  to  determine  1'0  at  the  start  of  a period.  The  fact  that  in  general  the 
response  to  a measurement  at  the  start  of  an  interval  does  not  occur  immediately 
constitutes  system  delay  time.  System  instability  will  occur  if  the  effect  of  this 
delay  is  not  counteracted. 
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The  wave  form  shown  in  Fig.  4 may  be  decomposed  into  the  sum  of  a square 
wave  and  a pulse  train  (see  Figs.  3.  3-2a  and  3.  3-3b). 


K T/2-H 


n/2-T0! 


Fig. 


3.3-3  Equivalent  Representation  ior  a Binary  Pulse  Pattern 


The  dc  value  of  the  wave  in  Fig.  3.3-3b  constitutes  the  essential  signal  feedback. 

Let  AT  =-  - T . The  digitizer  converts  the  analog  signal  w to  AT  by  the 
2 0 

following  relationship. 


AT  = -MIN(|w|/wmax, l)SIGN(w)  T/2 


(3-3) 


where  wMAX  is  the  analog  signal  corresponding  .0  tall  level  operation  (i.e. 
AT  = T/2). 


. i a and  the  torquer  constant  is  C (N/V), 

If  the  maximum  acceleration  is  ana  me  h v 


then: 


w — A M/C 
max  max  v 


(3-4) 


where  M has  been  previously  defined  as  pendulous  mass. 

The  wave  forms  of  Figs.  3.3-3a  and  3.3-3b  can  be  expanded  in  a Fourier 
series  of  period  T. 
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1'he  square  wave  is  given  by: 


-4/n  Z sin  (2m!  t/T)  /n 
r'odd 


l»-a) 


I lie  pulse  train  of  width,  T,  is  Riven  by: 


2 A T/T 


2/ „ I (-l)n+1  s in ( 2n r AT/T ) cos ( 2n” t/T ) /n 
[ n=l 


(3-6) 


The  total  signal  driving  the  torquer  is  then: 


2 AT /T  - P(T)  = 2 AT/T  - 4/n  Z sin  (2nnt/T) /n 

nodd 

- 2/n  Z (-l)n+1  sin  ( 2nn AT/T ) cos ( 2nn t/T ) /n 
n=  1 


(3-7) 


As  the  dc  term,  r.jy1  , is  proportional  to  w (see  Kq.  3-3),  the  periodic  functions 
of  T constitute  system  noise  of  known  period. 


With  no  input  (x  = 0)  the  signal  seen  by  the  digitizer  is  the  signal  of  llq.  3-7 
modified  by  its  propagation  through  the  dynamics  introduced  by  the  accelerometer, 
amplifier,  compensation,  etc.  The  steady  state  w at  any  sampling  time,  NT.  is 
given  by: 


= [ 2 AT/T  - P] A 


where: 


A 

u) 

0 


= (Mr/K)  c'oK0AMAX 
= 2 r/T 
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4>  - phase  shift  for  the  harmonic 

G(njQ>=  attenuation  of  the  n n harmonic 
P 1 Value  of  periodic  terms  at  sampling  time 

Substituting  for  AT  from  Eqs.  3-3  and  3-7  yields: 

w(NT)  = -Aim  + AC  / (A  M) 

V ITlaX 


MrC  K A P/K 
o o max 

1 + (rC  K C )/K 
o o v 


MA_  P/C 
max  v 


The  dc  value  of  wave  form,  w,  is: 


2AAT/T  =w(NT)  A/wmax=  A P 


= rC  K MA  P/(C  K) 
o o max  v 


(3-9) 


(3-10) 


Figure  3.3-4  shows  a typical  wave  form,  w,  for  zero  input  to  the  accelerometer. 


Fig.  3.3-4  Digitizer  Input  Waveform 

Note  that  from  Eqs.  3-9  and  3-10  the  ratio  of  w(NT), the  sampled  signal,  to 
the  dc  value  (essentially  the  maximum  amplitude  of  the  periodic  wave)  is: 


rC  K /K 
o o 

This  quantity  will  be  discussed  later  in  the  evaluations  of  the  four  digital 
i mple  ment  ati  ons . 
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I lie  compensation  network  shown  in  1 itl.  .1.2-2  has  the  following  transfer 
function: 

[T22s2  *■  aT2s  + l]tT42s2  + cT4s  + 1] 

[T32s2  + bT3s  + 1] 


1 lie  sampling  technique  (measuring  at  NT  and  respondin'  at  ( N * T ) introduces 
a phase  lug  equal  to  the  following: 


t>  (>  A _) 

where  a.  is  tile  frequency  of  the  signal  being  measured.  \s  the  mass  introduces 
180°  of  lag  shortly  after  its  5 second  breakperiod,  it  is  necessary  to  compensate 
for  this  delay  to  prevent  instability. 

1'he  lead  network  chosen  is: 

T 2s2  + cT.s  + 1 
4 4 


where: 


c = •**/  1 .8 

This  network  compensates  for  the  phase  lag  introduced  by  periods  greater  than  0.2 
second. 


The  remainder  of  the  compensation, 

[T22s2  +aT2s  + 1] 
[T32s2  +bT3s  + 1] 

where: 

T2  = 6/2  it,  a = 0.2 

T 3 = 12/2  n,  b = 2.0 
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was  introduced  to  attenuate  high  frequencies  (>2r/T),  thus  allowing  an  increase  in 
loop  gain. 

A plot  of  magnitude  and  phase  vs  frequency  for  open  loop  and  closed  loop 

response  is  shown  in  Figs.  3.3-5a,  b,  c and  d.  Linearity  was  found  to  extend  over 

09  -82  “3 

signal  levels  ranging  from  10  m/s  to  10  m/s  to  four  places.  Between  10 

“3 

and  5 x 10  accuracy  tapered  off  to  two  places  (see  Fig.  3.3-6). 


Amplitude  was  determined  to  be  independent  of  period  for  input  in  the  20  to 
40  second  range  (see  Figs.  3.3-7,  3.3-8  and  3.3-9). 

Figure  3,3-10  shows  the  effect  of  starting  out  in  saturation.  An  eight  second 
recovery  period  is  indicated. 

The  high  frequency  oscillation  (ringing)  occurring  during  the  first  eight  seconds 
in  Figs.  3.3-7  through  3.3-10  is  caused  by  the  transient  effect  of  starting  the  system. 
The  frequency  of  ringing  corresponds  to  the  frequency  at  the  -180  ordinate  in 
Fig.  3.3-5d. 


An  evaluation  of  the  binary  system  points  up  two  basic  weaknesses.  The  first 
which  is  encountered  in  all  sampled-data  systems  is  aliasing.  Aliasing  occurs  when 
the  signal  being  sampled  contains  frequency  components  greater  than  1/2  the  sampling 
rate.  This  condition  is  tolerable  if  only  a small  portion  of  information  lies  in  the 
proscribed  frequency  range. 


The  second  drawback  is  due  to  the  relatively  high  level  of  noise  in  the  form 
of  periodic  functions  existing  at  the  digitizer.  It  was  mentioned  earlier  that  the 
ratio  of  the  amplitude  of  the  periodic  signal  to  signal  for  no  input  is  rCQKQlK.  For 
the  parameters  considered  in  the  present  application, tCqKq/K.  ~ 10  . Due  to  the 
delay  inherent  in  the  operation  of  the  system.essentially  no  filtering  can  be  performed 
on  the  fundamental  of  the  periodic  functions  as  attenuation  introduces  intolerable 
phase  lag  (instability).  Although,  theoretically,  the  size  of  this  "noise"  creates  no 
problem  in  determining  the  signal,  hardware  considerations  preclude  the  accurate 
measurement  of  the  signal  in  the  presence  of  the  large, rapidly  changing, periodic 
function  (see  Chapter  5). 
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MACMTim:  (iiecib«l«) 


Fig.  3.  3- 5a  Closed  Loop  Magnitude  Response  for  Binary  Seismometer  System 


Fig.  3.  3- 5b  Closed  Loop  Phase  Response  for  Binary  Seismometer  System 
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TIME  <a) 

!•' ic.  3.3-10  Transient  Response  to  Saturating  Input 

3.4  Linear  System  with  Outboard  Analog-to- Digital  Conversion 

Outboard  analog-to-digital  conversion  implies  the  use  of  the  linear  system  in 
Fig.  3.1-1  with  a digitizer  attached  to  the  output  (see  Fig.  3.4-1). 


Fig.  3.4-1  Seismometer  System  with  Outboard  Digitizer 


For  simplicity,  the  feed  forward  components  have  been  lumped  into  Kj  and 
the  feedback  into  Kj.  A represents  a gain  of  1 for  the  digitizer.  Then: 


e_  = K,  A e,_/(l  + K,K,)  + nA 


(3-11) 


where: 


7j  = digitizer  noise 

e.  - input 
in 

eQ  = output 

Let: 

6A  represent  a change  in  the  digitizer  g;'  n 
6eQ  represent  the  resulting  change  in  the  output. 

T hen : 


6eo  = K:  ein  SA/(1  + K1K2)  + hSA 


(3-12) 


or  to  first  order: 

fie  /e  = fiA/A 
o o 


(3-13) 


i.e.,  the  fractional  change  in  output  is  equal  to  the  fractional  change  in  the  digitizer. 
3.5  Quasi-Linear  System  with  Inboard  Analog-to- Digital  Conversion 

The  A/ D converter  is  a quasi-iinear  instrument.  Measuring  infixed  increments 
or  quanta.  For  low  level  signals  (approximately  one  quantum)  the  output  is  not 
proportional  to  input.  A further  nonlinear  restriction  is  introduced  by  the  inability 
of  the  converter  to  track  rapidly  changing  signals.  This  constraint  may  be 
incorporated  into  an  analysis  by  considering  the  converter  a rate  limiting  device. 

The  above  analysis  applies  also  to  the  outboard  operation.  However,  there 
are  advantages  to  placing  the  digitizer  inside  the  loop,  advantages  which  usually 
result  from  feedback  networks.Figure  3.  5-1  shows  a configuration  with  the  digitizer 
inboard. 
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I'ijr.  :i.5-j  Srismomt'trr  System  with  Inboard  Digitizer 


lilt-  relation  between  output  and  input  is  given  by: 


o 

o 


K . A e.  / ( 1 + K.K.A)  + ’ A/ ( 1 + K.K-A) 
1 in  1 l 1 * 


where: 

A(nominal)  = 1 

It  is  apparent  from  Lq.  3- Id  that  noise  has  already  been  reduced  by  the 
(1  t K,K0A). 

Let: 

6 A = variation  in  digitizer  gain 
6e()  = effect  on  the  output. 

Then: 

? 2 

6e  = K,  e iA/(l  + K,K  AT  + h6A/(l  + Kj^A) 
o I in  i ^ 

To  first  order,  the  variation  in  output  is  given  by: 

6eo/eo  = (6A/A)/U  + kiX2A) 


(3-14) 


factor, 


(3-15 


(3-16 


Hence,  the  effect  or  the  percentage  variation  in  A is  reduced  b'  the  loop  gain. 


Computer  simulations  were  made  for  the  inboard  system  with  the  digitizer 
placed  in  the  high  gain  (.02-1)  Geotech  loop.  Step  inputs  were  used  rang.ng  in 
magnitude  from  4.0  x 10-10  m/s2  to  1.0  x 10’7  m/s2.  This  interval  was  chosen  to 
determine  the  effects  of  nonlinearities  on  system  response.  Overall  system  gain 
corresponded  to  the  Geotech  value  of  1.0x10  . FiveV/s  rate  limiting  in  the  digitized 
measurement  was  used  n all  the  simulations. 

Figure  3.5-2  shows  the  response  to  a step  input  of  4.0  x 10  ^ m/s.  Although 
the  effect  of  the  10  p.V  quantum  measuring  increment  is  evident,  the  output  has  the 
appearance  of  the  step  decay  appropriate  to  an  equivalent  linear  system. 


Gain:  105  V/(m/s2) 
Input:  4.  0 x 10"  m/s 
Quantum:  10.0  x 10  ® V 
Rate  Limit:  5.  0 V/s 


TIME  (s) 


Fig.  3.  5-2  Transient  Response  of  Inboard  Digitizer  to 
Low  Level  Step  Input 


with 


I igure  3.5-3  displays  the  response  to  a step  input  of  1.0  x 10  m s 
a quantum  of  40 /xV.  Again  the  figure  has  the  appearance  of  a step  decay. 

Figure  3.5-4  presents  the  res|>onse  of  the  system  to  a large  (relative  to  the 
40  /iV  quantum)  input  of  1.0  x 1 0*  ‘ m/s.  The  output  response  appears  linear,  as 
ttie  maximum  amplitude  is  several  orders  of  magnitude  greater  than  the  quantum 
of  the  digitizer. 

1'he  a cce'e ration  equivalent  of  a Rayleigh  displacement  wave  is  shown  ini*  igure 
3.5-5.  i lie  curve  was  derived  l>v  taking  second  differences  of  the  displacement 
data  and  using  linear  fits  for  continuity.  This  acceleration  equivalent  of  Rayleigh 
noise  generated  the  input  to  the  system. 

Figure  3.5-6  is  the  output  for  the  Rayleigh  acceleration  input.  A comparison 
of  the  two  figures  shows  good  agreement.  The  effect  of  the  rate  limiter  is  indicated 
by  tile  smoothing  of  large  rapid  changes  in  signal  ( 5 \ Is).  \ comparison  ot  the 
two  signals  was  made  at  the  acceleration  level  r ither  than  at  the  displacement  level 
since  this  approach  provided  a better  indication  of  system  performance.  Uouble 
integration  tends  to  smooth  out  disparities  in  comparisons. 
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Fig.  3.5-4  Transient  Response  of  Inboard  Digitizer 
to  High  Level  Step  Input 


} 
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ELECTRONICS  DESIGN  CONCEPT 

Within  the  framework  of  the  theoretical  analysis  presented  in  Chapter  3,  a 
basic  electronics  approach  was  chosen  to  meet  thf  desired  objectives.  A simplified 
block  diagram  of  this  approach  is  shown  in  Fig.  4.0-1.  With  the  A/D  converter 
removed  the  remaining  blocks  could  be  configured  as  in  the  current  Geotech  system. 


Fig.  4.0-1  Block  Diagram  of  Digital  Seismometer 

Control  Loop 


The  theoretical  results  indicated  that  a ternary  pulse-torque  loop  would  be 
unsatisfactory  and  further  loop  implementation  was  terminated. 

Since  both  the  pulse  width  modulated  loop  and  the  analog  loop  (A/D  converter 
within  the  loop)  required  equivalent  A/ 13  blocks  (through  the  pulsed  loop  required 
more  rigorous  constraints)  a study  was  made  of  the  A/D  portion  of  the  loop  in 
generalized  form. 

An  18  bit  A/D  converter  was  chosen  for  the  design.  The  full  range  ideally 
would  have  been  21  bits  to  correspond  to  the  seismometer  operating  goals  listed  in 
Chapter  2.  However , it  was  advisable  to  keep  the  resultant  hardware  from  becoming 
overly  complex  and  costly.  In  addition,  the  timing  requirements  of  the  pulse  width 
modulated  loop  (0.1  second  sample  period)  required  excessively  high  clock  frequencies 
in  the  A/D  converter  if  more  than  an  18  bit  A/D  converter  were  considered. 
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An  evaluation  o;'  the  noise  introduced  bv  the  A/I)  circuitry  is  presented  in 
Chapter  6.  This  noise  level  is  comparable  to  that  introduced  by  the  analog  |>ortion 
of  the  loop,  and  consistent  with  the  design  goal  of  having  the  instrument  noi<-e  .it  20 
db  below  the  lowest  anticipated  signals  to  be  measured. 
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CHART LR  5 

SLISMOMLTLR  DIGITIZING  LLLCTRONICS  ANALYSIS 

5.1  Introduction 

The  object  of  this  chapter  is  to  analyze  the  errors  and  problems  associated 
with  the  analog-to-digital  converter  in  a generalized  form.  The  analysis  splits  the 
A/D  converter  into  four  sections: 

(1)  \oltage  Reference  and  Scaler 

(2)  Ladder  and  Switch 

(3)  Amplifiers,  Sample  and  Hold,  and  Comparator 

(4)  Logic  and  Timing 

A representative  design  for  each  section  will  be  chosen  and  analyzed.  The 
errors  will  be  broken  into  static  errors  which  can  be  adjusted  to  zero  during  test 
and  compensatable  errors  such  as  temperature  coefficients  and  nonlinearities. 

5.2  Voltage  Reference  and  Scaler 

5.2.1  Voltage  Reference 

The  less  stringent  scale  factor  accuracy  requirements  for  this  A/D  converter 
relative  to  that  in  a gyro  loop,  make  it  possiole  to  amplify  the  basic  fi  V 1’VR  type 
reference  used  in  the  CSDL  SIRU  system  to  about  10  V,  using  an  operational  amplifier. 
This  will  introduce  some  static  errors  which,  however,  should  be  stable.  This  will 
also  allow  relatively  high  loads  on  one  reference,  since  the  output  impedance  will 
be  essentially  zero.  A 10  V level  was  picked  as  a good  compromise  between  the 
expected  input  signal  range  to  the  A/l)  converter  and  the  probable  use  of  a ±15  V 
bus  for  the  amplifiers.  Voltage  references  other  than  1’V  R's  could  be  used.  However, 
in  this  prototype  stage  it  is  advisable  to  use  a better  reference  and  then  back-off 
after  gaining  test  experience. 

5.2.2  Scaler 

There  are  two  ways  to  implement  the  requirement  for  reference  scaling  between 
the  lower  and  upper  order  bits.  One  is  to  scale  the  reference  as  applied  to  the 
resistance  ladders.  This  would  involve  a simple  resistor  divider.  The  second  is 
to  use  the  same  reference  voltage  for  both  ladder  sections,  but  scale  their  outputs 


k 
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by  the  gain  and  summing  resistors  in  the  comparator  section.  I'he  second  method 
is  more  desirable  since  it  allows  for  finer  adjustment  with  higher  value  resistors 
and  also  avoids  output  impedance  problems  associated  with  a reference  divider.  It 
also  places  all  the  critical  adjustments  required  for  the  A / 1)  converter  in  the  same 
location,  the  input  stage.  Thus  for  this  analysis  a single  value  reference  will  be 
assumed  and  scaling  will  be  done  in  the  comparator  stage. 

I'he  basic  implementation  is  shown  in  Fig.  5.2-1.  1 vpical  components  would 

be:  Vishav  type  S102  resistors,  Harris  HA  2904  operational  amplifier,  load 

compensated  SIKU  PVR  and  a good  grade  mylar  capacitor.  For  these  components: 

K^/R  - 1.6  (to  obtain  10  volts  output) 


e. 


Fig.  5.2-1  Voltage  Reference  Circuit 


Let: 

Rj  = 10  kO 

R2  = 16.0  k n 

Select  a value  of  C to  give  approximately  a 100  Hz  low  pass  filter.  This  will 

remove  some  of  the  noise  inherent  in  the  resistors,  amplifier  and  PVR.  The  required 

value  of  C is  given  by: 


C 


1 

2 irf  R 


.lgF 


(5-1) 
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5.2.3  Summary  of  Voltage  and  Scaler  Errors 


I 

I 


The  errors  for  the  voltage  reference  and  scaler  are: 
Static 

Tolerance  1.0  l^o 

Offset  57.1  /Ar 

Compensatable  ±20  pV  l°C 

See  Appendix  A for  a more  detailed  discussion. 


5.3  Ladder  and  Switch 

Ladder  networks  of  the  required  accuracy  are  available  only  to  12  bits.  Ihus 
the  ladders  and  switches  will  be  broken  into  two  sections  (see  Fig.  5.3-1).  I he 
main  problem  with  this  technique  is  to  maintain  a monotonic  summation  of  the  two 

ladder  outputs. 


i 


I 


I 


Fig.  5.  3-1  Ladder  Networks 


Figure  5.3-2  shows  the  outputs  of  the  two  ladders  graplm\dl> . The  ladder 
networks  are  divided  into  a least  significant  hit  12  (LSI!)  ladder  and  a (i  MSI!  ladder 
to  allow  adjusting  for  monotonicity  and  linearity  more  easily.  The  (i  hit  ladder  lias 
a small  resistance  subtracted  from  each  leg  to  allow  adjusting  for  the  "on"  resistance 
of  the  FET  switch  to  insure  monotonic  behavior. 


Fig.  5.3-2  Graphical  Representation  of  the  Ladder  Outputs 


Aj  and  are  gains  in  the  summation  amplifier  and  comparator  that  provide 
linearity  adjustment.  For  the  ladder  analysis, Aj  and  A^  are  assumed  to  he  unity. 
Gain  effects  will  be  discussed  in  the  next  section.  Overall  errors  will  depend 
primarily  on  the  initial  calibration  and  are  hard  to  estimate.  However,  the  calculations 
will  show  the  approximate  values  involved. 
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5.3.1  Monotonic  Requirements 

The  requirements  for  monotonicity  of  the  ladder  networks  can  be  examined 
as  follows: 

a)  6 MSB:  ASSUME  R (Fig.  5.3-1)  equals  10  kQ  and  Tj  equals  100  tK 

Assume  the  use  of  a Siliconix  DG190A  switch  with  a maximum  "on" 
resistance  of  30  Q ; also  that  RT  can  be  trimmed  to  within  0.1  Q over  a 
range  from  70  to  100  Q . 

Under  these  conditions  the  maximum  error  in  the  height  of  V 
will  be: 


1 2R_ 

10V  '2  “ 4R-r,+Rf_+ 


1 ' *'T  ^S-on 


10 


1 20,000 

[2  ~ 40,000  ± 0.1 


± 12 . 5pV 


The  height  of  one  step  in  the  12  LSB  will  be: 

HUS-  = 38.1viV 
218 

Thus  under  the  worst  case  of  the  highest  order  bit.  monotonicity  should  be 

attainable.  The  above  adjustment  also  compensates  for  the  initial  matching  tolerances 

of  R in  the  ladder  (typically  .01  to  .1%). 

b)  12 LSB:  ASSUME  R equals  10  kQ  and  r2  equals  30  fl.  Assume  the  same 

switches  are  used  as  for  the  MSB  ladder  and  that  the  RDS_on  match  is 

±10%. 


Under  these  conditions  the  maximum  error  in  LSB  height  will  be: 

10  ,1  20,000  s 

64  [2  ~ 40,000  ± 3 
± 5 . 8yV 

Thus  the  lower  12  bits  will  be  monotonic  without  further  adjustment. 


10  ,1 2R ) 

26  2 4R_r2+RDS-on 
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5.3.2  Temperature  Coefficients 

/■.••  • ni  (TCi  tracking  within  the  ladders  is  1 ppm/  C. 

•I'vnif-al  temperature  coefficient  irdLNl"i’ 

lyP  P . Tr  of  ,n  U\'I°C  on  the  MSB  of  the  ladder.  Because 

This  will  result  in  a worst  case  1C  of  10  w Co  (t50  ppm/*cl  has 

the  ratiometric  technique  is  used,  absolute  Thus  the  TC  of  the 

negligible  effect.  Switch  li,B.on  TC  is  typically  0.2  <1/  C. 

MSB  is: 


10  ,1 

64  ' 2 


20 , 000  t _ ♦ 5 . 8y  V 

a a a r o > 


40,000  + . 2 


• n nimhi  in%  therefore  monotonic 

The  switches  will  track  each  other  typically  to  • only  ,he 

behavior  will  be  assured  over  the  temperature  range.  • 

behavior  win  t h between  ladders  will  also  cause  a small 

overall  scaling.  2 ny  - resistors  should  introduce  no 

nonlinearity  in  the  overall  slope.  1 he  1C  of  the  trim 

significant  errors. 

5.3.3  Tolerance 

Tolerance  in  the  initial  value  of  11  in  Ihe  ladders  will  cause  only  second  order 
effects  since  the  R - 2R  technique  is  ratiometric. 

5.3.4  Leakage 

Typical  switch  leakage  a,  50°C  is  about  .4  nA  per  switch  The  ladder  never 
senses  these  leakages  because  the  SCOT  switches  shunt  them  to  REF 


(see  Fig.  5.3-3). 


DATA  6/r*  1 


DATA  B/T-0 


Fig.  5.3-3  Switch  Leakage  Paths 
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5.3.5  Switch  Transient 


The  worst  case  switch  transient  occurs  when  the  MSB  is  switched.  The 
magnitude  of  this  "spike"  is  determined  primarily  by  the  RC  product  of  the  ladder 
R and  the  switch  capacitance  (see  Fig.  5.3-4).  The  worst  case  transient  actually 
depends  on  the  initial  charge  states  of  the  various  capacitors,  but  for  this  approximate 
analysis  the  initial  conditions  are  assumed  zero.  The  time  constant,  i , can  be 
computed  by: 


T = + C + (L  ) X 4R 

T '^D  off  S on  D on' 

t = 19pF  x 40  kfi 

t = .76  us 


(5-2) 


Therefore  the  time  for  the  output  to  settle  to  less  than  1 LSB  error  (38  uV)  is: 


t 


t In 


V 

r REF 
LSB 


t = 9 US 


(5-3) 


This  time  would  limit  the  frequency  of  the  A/D  clock  to  about  110  kHz.  To  speed 
up  the  clock,  either  R of  the  ladder  must  be  lowered,  which  would  require  a more 
precise  adjustment  of  RT,  or  else  a shunt  should  be  placed  from  the  switch  common 
to  ground.  The  latter  method  is  more  practical. 


Fig.  5.3-4  Portion  of  Ladder  Network  Determining 
Switch  Transient  Response 
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Assuming  a 10  kQ  shunt  (l<s  in  P ig.  5.3-4),  the  time  required  to  obtain  less 
than  1 Lb B error  is  about  2.3  *r.  Therefore,  if  the  comparator  is  strobed  at  some 
time  greater  than  2.3  us  after  the  switch  command,  the  transient  would  introduce 
less  than  1 LSB  error  (Fig.  5.3-5). 


Fig.  5.3-5  Output  of  Ladder  Step 


5.3.6  Circuit  Delays  and  Maximum  Slew  Hate 


Circuit  delays  will  also  be  present.  These  will  consist  of: 


1 ) ladder  settling 

2)  switch  risetime 

3)  switch  deadtime 

4)  switch  drive  logic  risetime 


250  ns,  typical 
100  ns,  typical 
60  ns,  typical 
25  ns,  typical 


Total  delay 


435  ns,  typical 


The  total  delay  time  should  be  added  to  the  strobe  delay  time  giving  a total  strobe 
delay  of  about  2.6  us.  Assuming  the  strobe  time  is  step  centered  the  total  strobe 
delay  will  allow  step  times  of  about  5 us.  Therefore,  the  ladder  will  track  a maximum 
input  envelope  slew  rate  of  about  38  wV/5  us  = 1 Vis. 

5.3,7  Summary 

Because  of  the  adjustable  nature  of  the  ladder  and  switch  section  exact 
prediction  of  the  error  is  difficult.  However,  it  appears  that  the  circuit  and  components 
can  achieve  the  required  rate,  linearity  and  monotonicity.  It  should  be  noted  that 
the  necessity  of  ensuring  monotonic  behavior  sets  an  upper  bound  on  the  nonlinearity 
error  of  the  two  ladders.  The  remaining  nonlinearity  will  be  determined  by  the 
gain  errors  in  the  amplifiers  used  to  combine  the  two  outputs  and  the  temperature 
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coefficients  discussed  above. 


5.4  Amplifiers,  Sample  and  Hold,  and  Comparator 


5.4.1  Input  Amplifier 


The  input  amplifier  must  scale  the  input  signal,  E.n,  offset  it  to  allow  a single 
ended  output,  offset  any  static  errors  and  provide  a low  impedance  drive  for  the 
following  stages  (see  Fig.  5.4-1).  A high  frequency  rolloff  capacitor  is  also  shown 
to  allow  reduction  of  any  logic  generated  noise. 


MPUT  AMPLIFIER 


OUTPUT  TO 
( BINARY  LOOP) 
OR  COMPARATOR 
AMP 


Fig.  5.4-1  Input  Amplifier  Circuit 


Although  the  amount  of  offset  and  gain  necessary  in  this  stage  varies  with 
loop  configuration  it  will  not  exceed  unity  (-1)  gain,  and  me  offset  probably  will  not 
exceed  10  V.  Therefore  us  a worst  case,  unity  gain  and  10  V offset  will  be  considered. 
Offsets  of  less  than  10  V can  be  implemented  by  changing  the  (+)  input  gain  (adjusting 
R and  R.).  Precise  seismic  input  gain  is  adjusted  with  Rj  and  Rj.  The  high 
frequency  rolloff  occurs  at  frequencies  that  do  not  affect  the  stability  criteria  of 
the  loop,  and  it  need  not  be  computed,  since  its  exact  value  will  be  determined 

empirically. 

As  before,  an  HA2904  amplifier  and  Vishay  S102  resistors  are  assumed.  To 
prevent  loadingof  the  preceeding stages  inthe  seismic  amplifier,  Rj  isassumed  = 10 
kJl.  The  output  impedance,  ZQut.  may  be  computed  from  the  following  relationship: 
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7 = “ 

out  A 

1 + R_ 

(1  + ) (5‘4> 
R1 


where  A is  the  open  loop  gain.  For  the  HA2904,  A = 5 x 10  , typically,  and  /.q  = 200fl . 

1 herefore  /.  , a 0. 

out 

5. 4. 1.1  Slew  Hate 


The  maximum  slew  rate  of  any  expected  input  signal  is  96  \ sec.  (or  96aA  Iij s). 
The  11  \2904  slew  rate  is  specified  as  2.5  S 1/j s.  Therefore  no  error  is  expected. 

5. 4. 1.2  Summary  of  Input  Amplifier  Lrrors 

Assuming  unity  gain  and  10V  offset  the  errors  for  the  input  amplifier  are: 


Static  Errors 

Tolerance  .02% 

Offset  42  »V 


Compensatable  Krrors 
Signal  Gain 
Offset  Gain 
Amplifier  Offset 


i 1.5  ppm/°C 
± 1.5  [)j)m/°C 
± .8  »\'I°C 


As  in  the  reference  case,  the  static  errors  can  be  trimmed  out.  See  Appendix 
B for  a more  detailed  discussion. 


5.4.2  Sample  and  Hold 

In  the  case  of  a binary  force  rebalance  loop,  a significant  ripple  component 
will  be  present  on  the  input  signal  to  the  loop,  caused  by  the  10  Hz  current  pulses. 
Computer  simulation  has  shown  that  the  maximum  slope  of  the  ripple  component  at 
the  digitizing  time  is  96  V/s  (see  Fig.  3.3-4).  The  A/D  converter  cannot  follow 
this  fast  a signal,  even  in  a successive  approximator  configuration.  Thus  it  would 
be  necessary  to  sample  and  hold  the  input  signal  to  allow  an  A/D  conversion  if  the 
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binary  loop  were  to  be  used. 


If  a tracking  A/D  converter  is  used,  the  time  to  convert  depends  on  the  size 
of  the  change  from  sample  to  sample  with  an  envelope  slew  rate  of  4 V/s.  This 
change  can  be  .4  V,  which  is  equivalent  to  10,000  counts.  With  a clock  frequency  of 
200  kHz  a maximum  of  50  ms,  would  be  needed.  Clearly,  a tracking  configuration 
cannot  be  used  without  causing  large,  highly  variable  errors  and  loop  instability. 

Assuming  a successive  approximating  A/D  converter  is  used,  a maximum  of 
90//s  would  be  needed  to  convert  to  18  bits,  assuming  the  converter  resets  after 
each  conversion.  If  a tracking  successive  approximator  is  used  (i.e.,  the  A/D 
converter  remembers  the  previous  conversion  and  starts  from  thereat  each  sample), 
the  maximum  deviation  from  sample  to  sample  will  be  14  bits,  requiring  only  70  vs 
to  convert.  This  would  be  the  fastest  configuration  to  get  the  required  resolution 
on  a 4 V/s  envelope.  The  penalty  is  relatively  complicated  logic.  The  following 
analysis  will  attempt  to  design  the  sample  and  hold  and  to  analyze  its  errors  for 
this  configuration.  The  model  shown  in  Fig.  5.4-2  will  be  used. 


SAMPLE  MO  HOW  MODEL 


Fig.  5.4-2  Sample  and  Hold  Circuit 


Assume  the  turn-on  transients  have  decayed  io  a negligible  level  and  the  sample 
and  hold  is  sensing  only  a .4  V maximum  change  from  sample  to  sample  (every  0.1 
second).  Further  assume  the  sampled  input  value  is  10V.  In  order  for  the  error 
during  the  conversion  time,  rplL/Q « to  '3e  less  than  1 bit  or  38  ^V  (see  Fig.  5.4-3), 
the  RC  time  constant  required  is  computed  by: 


error | 


vln  - * RinCS> 


(5-5) 
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or 


R.  Cc  = 
in  S 


error. 

In  (1  - -y ) 

in 


Therefore: 


R Cc  ~ 18.4  seconds, 
in  a 


Assume  we  meet  this  requirement.  The  decay  can  now  be  calculated  by: 


decay  = V^n  d 


e 


t 


Rin 


(5-6) 


where: 


t - , 1 second 

V.  = 10  V 

m 

Therefore: 


decay  * . 054  V 

Thus  if  the  maximum  seismic  signal  change  from  sample  to  sample  is  .4  V, 
the  maximum  the  sample  and  hold  will  sense  is  .454  V.  The  sampling  aperture  is 
selected  so  that  less  than  1 bit  of  error  will  be  introduced  during  the  sample  time. 
Since  the  maximum  ripple  slew  rate  at  the  sample  time  is  96  V/s,  to  sample  with 
less  than  1 bit  of  error  the  sample  time,  rs>  would  be: 


38pV 

TS  - 96  V/s 


396  ns 


Let  t = 500  ns.  This  will  allow  for  the  rise  and  fall  times  of  the  switch. 
Using  the  computed  input  change  to  the  sample  and  hold  between  samples  (.454V), 
the  value  of  Cs  to  allow  charging  to  within  38  *V  (1  bit)  of  the  correct  value  in 
can  be  computed  (see  Fig.  5.4-4).  The  minimum  resistance  through  which  Cg  can 
be  charged  is  RDS_on.  Assuming  use  of  a DG190A,  RDS_on  is  30  Q.  Thus  Cg  can 

be  computed  by: 


| error  | = Win  e 


t 

RDS-onCS 


(5-7) 
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where: 


1 


M^n  = * 454V 


T 

S error 

Cs  rDS-ON1N  av 

Cg  - 1775  pF 


SAMPLED 

VALUE 


Fig.  5.4-3  Signal  Decay  During  A/D  Conversion  Time 


Fig.  5.4-4  Sampling  Error 


R.  is  then  computed  by: 
in 

R.  C*.  = 18.4 
in  b 

“in  * '°'°a 

This  high  value  Tor  R,„  is  impractical,  especially  for  a fast  enough  1.  C.  ampl.fier 
to  allow  a correct  comparison  to  be  made. 

It  should  be  noted  that  this  analysis  is  somewhat  optimistic  in  assuming  the 
ripple  component  slew  rate  does  not  change.  In  reality  this  will  not  be  the  case. 
The  96  V/s  value  was  generated  from  a simulation  assuming  i0%  duty  cycle  of  the 
binary  input  pulse.  For  duty  cycles  corresponding  to  large  dynamic  excursions. 
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the  slope  will  change  radically.  Furthermore,  effects  of  leakage  current  and  offset 
associated  with  such  high  values  of  H-n  would  also  introduce  large  errors.  Input 
offset  currents  of  the  comparator  alone  would  introduce  temperature  coefficient 
errors  of  greater  than  .5  V/°C  (or  > 13  bits /°C).  ttven  if  the  10U'  3 value  for  H.n 
could  be  achieved,  a large  nonlinearity  would  still  result,  thus  negating  any  advantages 
achieved  by  the  relatively  complicated  tracking  successive  approximation  technique. 


If  other  sampling  techniques  are  used,  such  as  peak  detectors,  the  loop  errors 
would  still  be  high  since  leakages  would  still  dominate  because  of  the  high  impedances 
necessary  to  obtain  the  accuracy  during  the  conversion  time.  A peak  detector  would 
also  force  the  A/I)  conversion  to  be  performed  after  the  occurrence  of  the  peak 
signal.  This  would  require  a conversion  to  occur  during  the  succeeding  force  rebalance 
pulse,  thus  restricting  the  dynamic  range. 

In  conclusion  the  excessive  speed  and  accuracy  problems  associated  with  the 
large  ripple  component  make  the  binary  force  rebalance  loop  extremely  difficult  to 
implement.  It  should  be  repeated  that  this  restrictive  sample  and  hold  design  is 
necessary  only  for  the  binary  rebalance  method.  Any  approach  not  having  the  large 
ripple  component  would  yield  a sample  and  hold  circuit  that  was  readilj  achievable 
(if  it  were  necessary  at  all). 

5.4.3  Comparator  Amplifier 

The  comparator  amplifier.  Figs,  5.4-5  and  5.4-6,  must  receive  the  offset 
adjusted  input  signal  (or  the  sample  and  hold  output)  and  compared  to  the  summation 
of  the  ladder  outputs  to  provide  a logic  command  indicating  either  the  equality  or 
the  larger  of  the  two  inputs.  It  also  must  have  a slight  amount  of  hysteresis  or 
dead  band  to  prevent  noise  induced  oscillation  about  its  final  value. 


SUMMING  AMR 


6 MSB  >2  LSB 

Fig.  5.4-5  Comparator  Amplifier:  Convert  and  Hold  Type 
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COMPARATOR 

SUMMING  AMP 


Fig.  5.4-6  Comparator  Amplifier:  Tracking  Type 

Two  configurations  are  possible  depending  on  whether  or  not  the  A/D  converter 
must  track.  Though  it  is  possible  to  combine  the  functions  of  the  summing  amplifier 
and  the  comparator  into  one  stage,  this  analysis  assumes  a separate  summing 

amplifier  to  allow  more  flexibility  of  adjustment.  For  both  amplifier  - comparator 

configurations  the  summing  amplifier  is  the  same. 


5.4.4  Summing  Amplifier 

Assume  a required  unity  gain  for  both  the  input  signal  and  the  combined  ladder 
input,  and  ' a scaling  of  the  12  LSB  input  to  the  6 MSB  by  a gain  of  1/2  or 

.015625.  Figv.  w j.4-7  then  shows  a generalized  differencing  amplifier.  Given  the 
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resistance  ratio  constraint  as: 


where: 


m 

l 

i = l 


k. 

1 


the  expression  for  v is  given  by: 


(5-8) 


vc  = 1 k'j  V'j  " Ki  Vi 

c j=l  J J 1=1 


(5-9) 


For  the  caseof  two  non-inverting  inputs,  one  signal  inverting  input  and  one  grounded 
inverting  input,  (shown  in  either  Fig.  5.4-5  or  5.4-6): 

vc  ■ ka  va  + kb  vb  - kin  vin  ' k0  ' 0 


or: 


v 

c 


R20 

( ) v . 

Rio  in 


(5-10) 


where: 


ka  + 


kb  = kin  + 


or 


The  use  of  the  grounded  inverting  input  resistor  is  necessary  to  allow  the  resistor 
ratio  constraint  to  be  soluble  for  unity  signal  and  ladder  gains. 

R 

For  unity  signal  gain,  ^ = 1 . The  corresponding  scaling  on  the  ladder  inputs 

is  therefore: 
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1 


R 

R 


£ 

a 


. 015625 


If: 


R 


20 


R 


C 


10  kft 


then: 

R20  = Ra  = Rc  " 10  kf2; 
640  kfl 

640  k 51 


5.4.4. 1 Output  Impedance  and  Slew  Rate 

As  in  the  case  of  the  input  amplifier,  ZQut  tt  0 and  the  slew  rate  is  high  enough 
to  contribute  negligible  errors. 

5. 4. 4. 2 Summing  Amplifier  Error  Summary.  The  errors  for  the  summing  amplifier 
are: 


Static 

Toleiance  .03% 

Offset  42  PV 


Compensable 

Signal  Gain 
Ladder  Gain 
Amplifier  Offset 

I 


±1.5  ppm/°C 
±1,5  ppm/°C 
± .8  ^ V/°C 


As  before  the  static  errors  can  be  trimmed  out.  See  Appendix  C for  a more 
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detailed  discussion. 


5.4.5  Comparator 

The  two  comparator  configurations  are  shown  in  Figs.  5.4-5  and  5.4-6.  Of 
the  two,  the  latter  will  have  the  higher  errors  because  of  the  need  for  two  comparators 
to  get  the  extra  logic  information.  Assuming  we  want  the  threshold  band  to  be  ± 1 
LiSB,  then  2Vj^  = 76  //V.  This  value  is  generated  by  the  dividers  off  the  supply. 
Let  Rjjj  = Kn3  = 1 M<?  . Then: 


2vd  " Ir 


RD2  + RD4 


D1  + RD2  + RD3  + RD4 


] x 30V 


2Vq  = 76tiV 


(5-11) 


Computing  yields: 


**02  2.4  fi 

Any  error  in  this  value  will  generate  only  a second  order  effect  in  the  comparator 
error  and  will  be  negligible.  The  deadband  will  also  be  easily  adjustable. 

Frrors  due  to  input  current  sources  will  also  be  negligible  since  the  output 
impedance  of  the  preceeding  stage  is  essentially  zero.  Errors  due  to  input  offset 
voltage,  although  high  in  magnitude,  can  be  nulled  out  by  proper  adjustment  of  R 
and  Rjj4,  or  by  offsetting  the  preceeding  stage.  It  is  assumed  that  LM  106 
comparators  will  be  used. 

5.4.5. 1 Summary  Comparator  Errors.  The  errors  in  the  comparators  are: 


Input  Voltage  Offset 

.5m V X 2 =1.0  mV 

Input  Offset  TC 

2 X 3 pV/°C  = 6 tN/°C 
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Response  Time 

The  response  time  of  the  LM  106  is 
typically  less  than  100  ns  for  2 mV 
overdrive  and  therefore  no  speed 
errors  will  be  introduced. 


5.4.6  Summary 

The  errors  introduced  by  the  amplifiers  and  comparators  all  are  well  within 
the  limits  allowable  to  design  an  accurate  18  bit  converter.  The  sample  and  hold 
errors,  nowever,  arenot.  Thus  any  binary  force  rebalance  design  should  be  avoided. 
Not  mentioned  in  any  of  the  above  calculations  are  power  supply  errors.  The  rejection 
ratios  of  all  the  components  used  are  high  enough  to  ensure  negligible  power  supply 
induced  errors  for  any  reasonable  power  supply  design. 


5.5  Logic  and  Timing 


Logic  and  Timing  circuits  for  this  application  do  not  present  any  difficulty. 
The  only  critical  area  in  the  design  is  in  the  speed  of  the  18  bit  counter  that  drives 
the  ladder.  For  any  tracking  approach,  low  power  TTL  or  CMOS  would  probably 
have  sufficient  speed.  If  a successive  approximator  is  used.the  counter  would  have 
higher  speed  requirements,  forcing  the  use  of  standard  power  TTL. 


For  the  binary  force  rebalance  approacn,  the  logic  would  be  more  complicated 
and  would  require  a high  frequency.quantizing  clock  at  about  2.6MHz  which  would 
have  to  be  crystal  controlled.  Other  approaches  do  not  require  frequencies  higher 
thanabout  200  kHz  and  the  frequency  stability  of  the  oscillator  used  would  not  affect 
the  accuracy  of  conversion.  Another  complexity  which  the  binary  technique  imposes 
is  the  need  for  a second  18  bit  counter  to  control  the  force  pulse. 


5.6  Summary 

Chapter  5 has  detailed  the  analog  error  sources  introduced  in  a generalized 
model  of  an  18  bit  A/D  loop  for  the  seismic  system.  Total  error  depends  on 
configuration;  however,  all  the  errors  seem  reasonable  as  long  as  a design  is  used 

which  does  not  introduce  large  ripple  components  on  the  analog  signal  to  be  quantized. 

The  binary  force  rebalance  scheme  is  rot  recommended  due  to  the  prohibitive  sample 

and  hold  design  requirement. 
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CHAPTER  « 


NOISE  ANALYSIS  Ol'  AN  A LOO- TO-  DIGIT  A L CONVERT  Ell 
6.1  Introduction 

The  loop  used  for  this  analysis  is  the  analog-to-digital  converter  (A/D) 
(discussed  in  Chapter  5).  It  should  be  noted  that  noise  estimates  are  highly  dependent 
on  the  circuit  configuration;  therefore,  this  estimate  is  valid  only  for  the  loop  topology 
selected,  shown  in  Figure  6.1-1, 


n. 


Fig.  6.1-1  Analog-to-Digital  Converter 
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ti.2  Definition  of  Noise  Sources 
6.2.1  Resistor  Noise 

Nyquist  noise  generated  by  the  resistors  is: 

e = / 4kTBR 
n 


or 

e = 1.27  x 10’10  /"BR  rms 
n 

where  B is  the  effective  bandwidth  and  R is  the  resistance. 

The  highest  bandwidth  encountered  in  the  design  is  10  Hz  and  the  maximum 
effective  resistance  is  less  than  lOk^;  therefore,  the  maximum  expected  noise  from 
an^  resistor  is: 


e = 1.27  x 10"10  /105 
n 


e = . 04pV  rms 
n 


This  value  is  negligible  compared  to  other  sources.  Thus,  the  resistances  in  the 
circuit  of  Fig.  6.1-1  will  be  assumed  noiseless. 

6.2.2  Reference  and  Buffer  (VREp) 

Using  the  notation  on  Fig.  6.1-1  and  the  subscript  n to  indicate  the  noise 
component: 

VP  = (VPVR  + VPVRn) 


or: 


where: 


VREF  = - zr-  Vn  ‘ <1  + 5^)  A 


R1  P 
Z2  Z2 

VREF  = R^  VPVR  + R^  VPVRn 


R_ 


R1  in 


(1  + R?  Aln 


= r i I 

lr2c2s+1j 


(6-1) 
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6.2.3  Input  Amplifier  (v„J 


IN' 

R r7  r7 

VIN  = R^  VREF  ' R^  EIN  + ( 1 + R^'  A2n 


Substituting  Eq.  6-1  into  the  above  yields: 


IN 


R5Z2  R5Z2 

RxR4  VPVR  RxR4  PVRn 


R r Z2 

RT  (1  + r7>  Aln 
4 1 


R7  R7 

+ (1  + 


2n 


6.2.4  Summing  Amplifier  (vc> 


v = — V 
a 26  REF 


M ,, 

Vb  ~ ^T2  VREF 


R12  R13  R13  Rx? 

vc  = R^  VIN  + -R^  Va  + -R^  Vb  + (1  + rT:]  A3n 


'10 


Reducing  these  equations  yields: 


R12  rR13  N R13  M 

VC  * R^  V!N  + L R9  26  „8  2 


r]  Vr 


+ (1  + 


12 


)A, 


R„  „T7J  VREF  T ’ R^'"3n 


6.2.5  Total  Noise  Equation 

Substituting  Eqs.  6-1  and  6-2  into  this  expression  yields: 


R12  R5Z2  „ 

Vc  " R1q  R,R4  PVR 


R12  R5  Z 2 
R10  R1  R4 


PVR. 


n 


R-ijRc  R12R7 

. 12  5 ,i  R j.  _l£ — !_  p 

*10*4  R1  ln  IN 


'lo1^ 


R12  r7  Z2  r 1 
' sn  IL  * ^ 2n  * R1  R 


R13  N 


'13 


M 


R 


13 


'10 

z,  r13  N 
+ _£  + 

Rx  L R9  2 6 

Z_  R,,  N 

- [1  + /]  [^—6  + 
R1  K9  2b 


9 2' 


R8  212 


M 


"“T?11  VpVRn 


R 


13 


M 

T2J 


R 


PVR 


12 


_ i,-]  A.  + d + R 

Rg  212  In  R10 


) A 


3n 


(6-2) 


(6-3) 
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As  can  be  seen  above,  several  of  the  noise  terms  are  introduced  in  phase 
with  each  other  and  can  be  summed  before  the  noise  components  are  root  sum 
squared.  Grouping  terms  yields: 


v = 


R12  R7 
*10  R6  IN 


2 r ,R13  N_  + R13  M ,)  _ V 

' ~ IT'  " 


**T[' 


R9  2 6 R8  212 


*10  V PVR 


Z2  r,R13  N . R13  M , _ R12_R5]  „ 

R7  [(R^6  + R10  V PVRn 


'9  2 


- U + R^>  [<R 


R13  N 


R13  M 


9 2 6 R8  212 


) - 


R12  R5i  a 
R10  *4  ln 


(6-4) 


R1 2 R7 

— (1  + vr-)  A- 

R10  R6  2n 


+ (1  + 


'12 

l10 


) A 


3n 


The  first  two  terms  in  Eq.  6-4  are  the  signal  terms,  while  the  remaining 
terms  are  noise.  We  can  now  calculate  the  various  output  components  of  the  noise 
due  to  the  various  noise  sources. 


6.3  Assumptions 

The  following  assumptions  are  made  before  evaluating  the  noise  defined  by 
Equation  6-4: 

1)  Neglect  resistor  noise 

2)  Neglect  ladder  noise 

3)  Assume  we  place  a 10  Hz  filter  on  vc  that  introduces  no  noise.  Thus 
the  comparator  sees  only  the  noise  components  below  10  Hz.  The  filter 
is  placed  at  this  point  in  the  system  because  it  is  the  most  effective 
position  to  reduce  the  noise  components  introduced  by  the  amplifiers. 
If  filters  are  used  on  the  feedback  elements  in  amplifiers  A2  and  A3, 
only  part  of  the  noise  introduced  by  A2  and  A3  goes  through  the  feedback 
filter.  The  remaining  passes  through  with  unity  gain  and  no  bandwidth 
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limiting.  In  addition,  several  of  the  noise  sources  would  no  longer  cancel 
because  of  phase  delays. 

4)  l’he  filter  on  the  l’VH,  Z2>  is  assumed  to  rolloff  at  10  Hz.  This  filter 
can  be  used  in  the  feedback  loop  of  A1  since  phase  delay  here  does  not 
alter  the  cancellation  of  noise  in  later  stages,  though  still  reducing  a 
part  of  the  noise  on  Vp^jrp  above  10  Hz. 

5)  The  resistor  values  and  components  determined  in  the  proceeding  sections 
* ill  be  used.  These  are  by  no  means  optimum. 

6.4  Noise  Computation 
6.4.1  1’VU  Noise  Component 

Noise  due  to  PVR  = 

h.  rA3  N + Ri3  M_ ) 

R1  R9  26  R8  212 

Measured  data  on  S1RU  type  l'VB's  gives  values  of  about  1.8  pV  rms  for  dOO 
Hz  bandwidth.  Over  a 10  Hz  bandwidth  1.0  #A’  rms  would  seem  a reasonable  estimate, 
since  there  is  some  1/f  noise  present. 

The  worst  case  condition  of  noise  due  to  the  PVR  occurs  when  full  mass  position 
offset  is  introduced  (R4  = R,.,  and  the  digital  control  is  supplying  N 8 M - 0). 
For  this  case. 


R,,  Rr 

_ Lr__ 2.1  v 
R10  V PVRr 


Noise  due  to  PVR  (worst  case): 


= 1.6x  10 

= 1,6  uV  rms  for  10  Hz  bandwidth. 
6.4.2  A1  Noise  Component 


Noise  due  to  Al: 


[( 


R13  N 


Z + R 


R13  M 


A 


In 
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Worst  case  noise  occurs  under  the  same  conditions  as  for  the  PVR,  where  N M 0 
and  R^« 

6.4.  3 IIA  2904  noise  is  specified  as  dominantly  voltage  noise  of  900  n^z  at  10Hz. 
therefore: 


Aln  = 2.84  pV  rms  A2n  = A3n) 


Thus: 

Noise  due  to  Aj  = (1  + 1.  6)  x 1 * 2.  85  pV  rms 
= 7.41  pV  rms 

(for  10  Hz  bandwidth) 


6.4.4  A2  Noise  Component 

R12 

Noise  due  to  A„  = (1  + a ) ^2n 

= 1 x 2 x 2.  85  pV 
= 5.  70  pV  rms 

(for  10  Hz  bandwidth) 


6.4.5  A?  Noise  Component 

Noise  due  to  A^ 


(1 


12 

10 


) A 


3n 


2 x 2.85  pV 

5.  70  uV  rms 

(for  10  Hz  bandwidth) 


6.  4.  6 Total  Noise  - Worst  Case 

The  total  noise  at  the  output  is  the  root  sum  square  of  the  above  sources 

Total  noise:  + 7.412  + 5.  702  + 5.  702 

Total  noise  = 11.1  uV  rms  for  10  Hz  bandwidth. 
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6. 5 Summary 


The  derived  value  of  11.1  pV  rms  is  the  worst  case.  In  a more  realistic 
case  only  a small  amount  of  offset  will  be  used.  For  that  case  the  noise  level  will 
be  dependent  on  the  digital  control  number  driving  the  ladder.  Figure  6.4-1  shows 
a graph  of  noise  at  the  comparator  versus  the  digital  number,  assuming  various 
mass  position  offsets. 

It  should  be  noted  that  this  noise  estimate  does  not  take  into  account  the  seismometer 
noise  or  the  effect  of  the  rate  at  which  the  A/ 1)  converter  output  is  sampled.  Roth 
of  these  will  alt  act  the  overall  system  noise. 


Pig.  6.4-1  Output  Noise  vs  Fraction  of  F'ull 
Scale  Ladder  Voltage 
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CHAPTER  7 


i 


CONCLUSIONS  AND  RECOMMENDATIONS 

» 

I. 

7.1  Conclusions 

| 

T \e  present  study  examined  three  separate  approaches  to  implementing  a 
feedback  loop  around  the  Geotech  36000  seismometer  while  providing  a digitized 
output.  The  first  approach,  using  a ternary  pulse  torque-to-balance  loop  was  shown 
analytically  to  result  in  a generated  noise  signal  (caused  by  the  pulse  torquing) 

which  was  input  dependent  and  could  not  be  removed  effectively  from  the  desired  I 

data  signal. 

The  second  approach,  using  a pulse  width  modulated  loop,  eliminated  the 
problem  encountered  in  the  ternary  implementation  by  constraining  the  pulse 
frequency  to  be  constant  and  far  enough  removed  from  the  desired  data  signal  so 
» as  to  permit  a stable  loop  with  good  linearity  and  acceptably  low  noise.  However, 

the  pulse  torquing  gave  rise  to  a large  ripple  signal  riding  on  the  output,  lhis 
ripple,  though  tractable  analytically,  was  of  such  magnitude  that  it  produced 
insurmountable  sampling  problems  in  the  electronics. 

i The  third  approach,  using  an  analog  loop  with  A/D  conversion,  proved  to  bo 

the  best  alternative.  Of  the  two  possible  implementations  of  such  a loop,  inboard 
i or  outboard  of  the  feedback  loop,  the  inboard  implementation  appears  to  be  much 

superior. 

The  prime  disadvantage  of  the  external  A/D  converter  is  that  its  errors  add 
directly  to  the  other  instrument  errors  and  are  not  easily  amenable  to  determination 

and  compensation.  ' 

Using  an  A/D  converter  within  the  loop,  as  discussed  in  Chapter  3,  retains 
the  simplicity  of  the  simple  add-on  external  approach.  But  most  important,  almost 
all  errors  resulting  from  the  A/D  are  reduced  by  the  gain  of  the  feedback  loop, 
allowing  for  a more  accurate  and  easily  fabricated  circuit. 

The  worst  case  noise  resulting  from  the  A / D conversion  was  found  to  be  slightly 
higher  than  that  generated  within  the  seismometer  (11  vV  rms  (0-10Hz)  vs  4 /A 

rms  (.02-1)).  The  overall  system  noise,  with  the  inboard  A/D.can  be  expected  to  be  ^ 

comparable  to  that  with  the  present  analog  system.  | 
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For  both  inboard  aid  outboard  implementations, a relatively  slow  A/D  can  be 
used  since  the  digitized  signal  is  a "mirror"  of  the  seismic  signal  and  has  the 
same  low  rate  of  change.  There  are  no  further  sampling  constraints  at  frequencies 
of  interest,  thereby  reducing  any  aliasing  problems.  In  short,  from  both  analytical 
and  electronics  standpoints  the  use  of  an  A/D  converter  inside  the  loop  appears  as 
an  attractive,  easily  implemented,  digital  solution. 

7. 2 Recommendations 

Ilased  on  the  results  of  the  foregoing  study  the  following  recommendations 
are  presented: 

1)  An  inboard  A/D  converter  should  be  interfaced  with  a Geotech  36000 
seismometer  and  the  resultant  performance  evaluated  at  a "seismically 
quiet"  site. 

2)  Comparison  should  be  made  of  the  seismometer  performance  for  a totally 
analog  configuration  (as  presently  used)  with  that  for  theabove  described 
inboard  A/D  configuration. 

3)  An  18  bit  A/D  converter  should  be  used.  A larger  range  will  reauire 
much  greater  cost  and  complexity  in  the  electronics  design.  Such 
measures  do  not  appear  justified  until  more  instrument  performance 
data  is  available. 
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APPENDIX  A 


VOLTAGE  REFERENCE  AND  SCALER  ERRORS 


V ditit  TOLERANCE: 
REF  


Resistors:  2 x .005% 

PVR:  1% 


T otal 


01% 

1% 


1.  01% 


VREF  OFFSET: 


HA2904 : 

voff  x (1  + 1.6) 


= 20  MV  x 2.6 


Total  I x 16000  xl.6=  2x10 
off 

T otal 


-10 


x 16000  x 1. 6 = 


52  yV 
5.  luV 
= 57.  luV 


V TC- 

V REF 


R^tg  tracking:  11.5  ppm/°C  x VpvR 


PVR:  tl  ppm / C 

HA2004 : 


AVoff  x (1  + 1.6)  = .4  Vrc  x 2.6 

AIoffX  16000  = 3 x 10"12  x 16000 
Total 


i 9.3UV/°C 
1 10.0  UV/°C 


t 1.0uV/°C 
♦ 0.05tV/°C 
20  uV/°C 
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APPENDIX  B 


INPUT  AMPLIFIER  ERRORS 


Signal  Gain  Tolerance: 

2 x .005% 

Offset  Gain  Tolerance: 

2 x .005% 

Total 


.01% 


.01% 

.02% 


Output  Offset: 


HA2904 

V x (1  + 1)  = 20  uV  r 2 

°ff  4 -10  4 

IQff  xlO  =2x10  1 x 10 

Total 


40 

2 

42  uV 


Signal  Gain  TC: 


KjHj  1. 5 ppm/°C  tracking 

Offset  Gain  TC: 

RgR^  1. 5 ppm/°C  tracking 

Amplifier  Offset  TC: 

HA2904 

AVQff  x (1  + 1)  = .4  pV /°C  x 2 
Aloff  x R4  = 3 x 10‘12A/°C  x 104 

Total 


±1.5  ppm/°C 


± 1 .5  ppm/°C 


± .8  pV/°C 
± .03  pV/°C 
± .8  pV/°C 
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APPENDIX  C 


SUMMING  AMPLIFIER  ERRORS 


Signal  Gain  Tolerance: 

3 x .005% 

Ladder  Gain  Tolerance: 


3 x .005% 

Output  Offset  (HA2904): 


off 

!off 


x (1  + 1)  = 20  yV  x 2 

-in  4 

x 10k  = 2 x 10  1U  x 10 
Total 


Signal  Gain  T.C.: 


.015% 


.015% 


40  uV 
2 UV 
42  uV 


R0'R20'R10 

TC  Tracking:  1.5  ppm/  C 
Ladder  Gain  TC: 

Ra'Rb'Rc 

TC  Tracking:  1.5ppm/°C 
Amp  Offset  TC  (HA2904): 

AV  ,,  x (1  + 1)  = .4uV/°C  x 2 

Oil  a 1 9 4 

AIoff  x 10  = 3 x 10  A/°C  x 10  ft  = 

Total 


±.8  pV/°  C 
+ .03  uV/°C 
±.8  uV/°C 
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